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Qi Abstract 
U 

We numerically investigate the impact of the General Theory of Relativity (GTR) 
on the orbital part of the satellite-to-satellite range p and range-rate p of the twin 
GRACE A/B spacecrafts through their post-Newtonian (PN) dynamical equations 
of motion integrated in an Earth-centered frame over a time span AP = 1 d. The 
present-day accuracies in measuring the GRACE biased range and range-rate are 
a p ~ 1 — 10 pm, Up ~ 0.1 — 1 pm s . The GTR range and range-rate effects turn out 
to be Ap = 80 pm and Ap = 0.012 pm s" 1 (1PN gravitomagnetic), and Ap = 6000 
pm and Ap = 10 pm. s _1 (1PN gravitoelectric) . It turns out that the range shifts 
Ap corresponding to the GTR-induced time delays At on the propagation of the 
electromagnetic waves linking the GRACE spacecrafts are either negligible (1PN 
gravitomagnetic) or smaller (1PN gravitoelectric) than the orbital effects by about 
1 order of magnitude over AP = 1 d. We also compute the dynamical range and 
^ 1 range-rate perturbations caused by the first six zonal harmonic coefficients J^,£ = 



2, 3, 4, 5, 6, 7 of the classical multipolar expansion of the geopotential to evaluate 
their aliasing impact on the relativistic effects. Conversely, we quantitatively, and 
preliminarily, assess the possible a-priori "imprinting" of GTR itself, not solved-for 
in all the GRACE-based Earth's gravity models produced so far, on the low degree 
zonals of the geopotential. The present sensitivity analysis can also be extended, in 
principle, to different orbital configurations in order to design a suitable dedicated 
mission able to accurately measure GTR. Moreover, it may be the starting point 
for more refined numerical investigations concerning the actual measurability of the 
relativistic effects involving, e.g., a simulation of full GRACE data, including GTR 
itself, and the consequent parameters' estimation. Finally, also other non-classical 
dynamical features of motion, caused by, e.g., modified models of gravity, may be 
considered in further studies. 
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1 Introduction 



Th e Gravity Rec o very and Climate Exp e rimen t (GRACE) mission (ITapley and Reigber . 



2001 




Taplev et al. 




2004 





Schmidt et al.l . 120081 ). jointly launched in March 



2002 by NASA and the German Space Agency (DLR) to map the terres- 
trial gravitational field with an unprecedented accuracy, consists of a0 tan- 
dem of two spacecrafts moving along low-altitude, nearly polar orbits (see 
Table [T] for their orbital parameters) continuously linked by an inter-satellite 
microwave K- band ranging (KBR) system accurate to better t han 10 /xm (bi- 
ased range p) (Reigber et al. . 20051 ) and 1 /im s _1 (range-rate) ( Reigber et al. 



20051 ; Kim et al.l . 1200 ll ). Investigations concerning a foll ow-on of the GRACE 



-fl- 

mission are being currently performed (IWei et al.l . 120091 ) ; by using an interfer- 
ometric laser ranging system it would be possible to reach an accurac y level 



of ~nm s 1 or better in measuring the range-rate ( iLoomis et al.l . 120061 ) . 



Although GRACE was not specifically designed to directly test the Gen- 
eral Theory of Relativity (GTR), which, indeed, was never solved-for in the 
several global gravity field solutions^ produced so far by different institu- 
tions from long data records from GRACE, the great accuracy in its KBR 
may, in fact, allow, in principle, to measure some consequences of GTR by 
expl oiting such direct accurate observables. Concerning the Lense-Thirring ef- 
fect (ILense and Thirringl . Il918l ). connected with the rotation of the source of 
the grav itational field (see eq. fl2| and e g. ([3]) below) , a similar idea was envis- 
aged by IChicone and Mashhoonl (120061 ). Thus, it is important to investigate 
the impact of the dynamical effects of GTR on both range and range-rate to 
preliminarily check if it falls within the present-or future-sensitivity domain 
of GRACE-type missions. It will be the subject of Section [2] in which both 
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3 See|http://icgem.gfz-potsdam.de/ICGEM/ on the WEB. 



2 



numerical and analytical calculations will be performed. To this aim, we stress 
that we will deal in depth with the effects of GTR on the range and range- 
rate of GRACE coming only from the orbital motions of the twin spacecrafts. 
Indeed, as shown in Section m the range perturbations corresponding to the 
GTR time delays affecting the propagation of the electromagnetic waves mu- 
tually linking GRACE A and GRACE B are not par ticularly important for 



our purposes. About the Lens e -Thirring time delay ( jKopeikinl . Il99ll . 11997 



Kopeikin and Mashhoonl . |2002| ; ICiufolini et al.l . 120031 ) . in Section 14.11 we will 



show that it is completely negligible for GRACE. The Shapiro time delay 
(jShapirol . Il964 ) is, instead, large enough to be detectable, but the magnitude 
of the corresponding range shift is smaller by about one order of magnitude 
than that due to the orbital motions (Section 14.21) . It should also be noticed 
that we do not aim to quantitatively assess the actual measurability of the 
post-Newtonian dynamical effects investigated. It is a different and important 
task which would deserve a dedicated work. Indeed, a fit of the initial condi- 
tions with quite a lot of other dynamical and empirical parameters to the real 
(or realistically simulated) observations would be needed in order to plausibly 
evaluate the level of removal of the effects of interest from the signatures. It 
is a nontrivial task which is beyond the scope of the present analysis. Con- 
cerning other, clock-related relativistic effects in GRACE, it must be recalled 
that dual-frequency carrier-phase Global Positioning System (GPS) receivers 
are flying on both satellites. They are used for precise orbit determination 
of both the GRACE A/B spacecrafts, and to time-tag the KB R system; the 



relatiy istic effects in the GRACE GPS data were examined by lLarson et al. 
fl2007l ). 



Some of the Earth's gravity field solutions retrieved from GRA CE data were 



used as background ref erence models in the LAGEO S-based tests f|Ciufo 



2009 



ini et al.. 



Iorio et al.l . l201ll ) of the Lense-Thirring effect (ILense and Thirringi . Il918l ) ; 

the even (£ = 2k, k = 1, 2, 3, ...) zonal (m = 0) harmonic coefficients Je of the 



multipolar expansion of the classical part of the terrestrial gravitational po- 
tential, estimated as solve-for parameters in the GRACE-based models, may 
retain an a-priori "imprinting" by GTR itself which, as already noted, has 
never been explicitly solved-for so far in the GRACE data processing. This 
aspect will be treated in Section [3J Section \5\ is devoted to the conclusions. 



Finally, let us remark that the approach followed here in the specific case of 
GTR can well be extended to other dynamical effects predicted, e.g., by mod- 
ified models of gravity. A first example is the investigation of the effects of a 



Yukaw a-like extra-force on the orbit of GRACE- A performed by lHaranas et al 



(120111). Such a perspective h as recently gained a new appeal after the proposal 
by lDvali and Vikmanl (120121 ) about a putative Earth-sourced fifth force as pos- 
sible explana tion of observed n eutrinic phenomenology in some Earth-based 
laboratories ( lAdam et al.l . l201lh . 
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2 General relativistic effects in the satellite-to-satellite range and 
range-rate 



The relativistic treatment of the near-Earth s atellite orbits, as far as the 
equations of motion are concerned, is outlined in iPetit and Luzuml ( 120101 . p. 
155) and references therein. 



2.1 Numerical calculation 



The numerical approach adopted is as follows. 



We simultaneously integrated the equations of motion of both the GRACE 
spacecrafts in a geocentric local inertial frame3 endowed with cartesian coor- 
dinates by including the GTR effects we are interested in as dynamical orbital 
perturbations of the Newtonian monopole. More specifically, the Parameter- 
ized Post-Newtonian (PPN) perturbing acceleration Appn to be added to the 
Newtonian monopole ANcwton = —GMf/r 2 in the equations of motion 



alt 2 



L Ncwton 



LPPN 



id^J, to order 0{c~ 2 ) (1PN), (ISoffei Il989l . p. 89, p.95), (IPetit and Luzuml . 

(2) 



20101 . p. 155) 



l PPN 



-23,-2 



x B 



g- 



In it (ISoffell . Il989l p. 89, p.95), (IPetit and Luzuml . l2010l p. 106), 



E n 



GM f 

c 2 r 3 |^ 



2(/3+ 7 )GM _ 2 ' 



r + 2(1 + 7) (r ■ v) v } 



(3) 



B - - fi±2) ^ 
n* — I 2 ) cr3 



3(S-r) 



where 7,/? are the usual PPN parameters (jWilll . Il993l ). c denotes the speed 
of light in vacuum, G is the Newtonian constant of gravitation, M and S are 
the mass and the proper angular momentum, respectively of the central body, 
and v is the velocity of the test particle moving at distance r from M; the unit 
vector f is directed from the central body to the test particle. In eq. ((2])-eq. ([3]) 
Eg is the so-called "gravitoelectric" , Schwarzschild-like fi eld, while B g is the 



"gravitomagnetic" one yielding, among other things, the iLense and Thirring 
(119181 ) orbital precessions. We worked in the GTR case, i.e. for 7 = (3 = 1. 



4 We neglected non -inertial effects caused in it by t idal forces due to external bodies 
of the solar system ( Brumbeg and Kopeikin . 19891 ). 

5 It is assumed that the de Sitter precession has been transformed away. 
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Then, we performed another integration without them. 



The time span of both the integrations was AP = 1 d, in order to keep con- 
tact with the actual data re duction procedures followed by the GRACE ana- 
lysts ( jFlechtner et al.Ll201lh . The method adopted, implemented with the soft- 
ware package MATHEMATICA, is the ExplicitRungeKutta one, with MaxSteps — > 
10 6 and MaxStepFraction — » 1/1000. The initial conditions, common to all 
the numerical integrations, are in Table [TJ The altitudes of the twin GRACE 

Table 1 

Keplerian orbital elements of the pair GRACE A/B correspond- 
ing to the state vectors, in cartesian coordinates, of the files 
GNVlB_2003-09-14_A_00 and GNVlB_2003-09-14_B_00 retrieved from 
|ftp://cddis.gsfc.nasa.gov/pub/slr/predicts/current/graceA_irvs_081202_0.gfz] and 
ftp:/ / cddis.gsfc.nasa.gov / pub / sir / predicts / current / graceB_irvs_081201_l .gfz. The 
epoch is 13 September 2003. See ftp:/ /podaac.jpl. nasa.gov/pub/grace/doc/Handbook_lB_vl. 3.pdf 
for the explanation of the GPS Navigation Data Format Record (GNV1B) format. 
The orbital elements are the semi-major axis a, the eccentricity e, the orbital 
inclination / to the Earth's equator, the longitude of the ascending node 0, the 
argument of pericenter lj, and the mean anomaly Ai. The Keplerian orbital periods 
p(Kep) ^ 27r v / o 3 /GM ffl of the GRACE pair are of the order of » 1.56 h= 0.065 d. 



s/c 


a (km) 


eo 


h (deg) 


n (deg) w (deg) 


M (deg) 


A 


6841.11877 


0.00272831 


89.9395 


-71.5742 119.916 


-179.997 


B 


6839.80210 


0.00298412 


89.8374 


-71.5081 118.082 


-179.997 



spacecrafts are about 500 km with respect to the Earth's surface; their orbits 
are almost circular and polar. 



The resulting numerically integrated trajectories were, then, used to com- 
pute the satellite-to-satellite range perturbation Ap as the difference among 
the perturbed and the unperturbed ranges. The range-rate perturbation Ap 
was straightforwardly computed by numerically differentiating Ap. 



2.2 Analytical calculation 



Our numerical approach was successfully tested by comparing its outcome 
for the Lense-Thirring effect to an analogous analytical calculation of the 
gravitomagentic range and range-rate signals. Here we outline it in detail. 

As a first step, it is required to work out the short-period, i.e. not averaged 
over one full orbital revolution, Lense-Thirring perturbations of all the six 
osculating Keplerian orbital elements a, e, J, Q, u, M. of a test particle. It can 
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1989 



1991 



be done in the framework of the standarcfjll R — T — N formalism ( Soffel , 
p. 90) a pplied to the perturbative Gauss equations (IJoos and Grafarend . 
p. 25), ( ISoffell . Il989l p. 90) for the variation of the orbital elements. B y using 
the R — T — N components of the gravitomagnetic force ( ISoffell . Il989l . p. 95), 
( IJoos and Grafarendl . Il99ll . p. 24) in eq. (j2J) , ^mputed in a frame with the z 



axis directed along S, it is possible to obtain^ 



Aa LT 



Ae 



LT 



AI 



LT 



AQ 



LT 



cos IAQ LT + Awlt 



AM 



LT 



2GS cos /(cos /—cos /o) 
c 2 na 3 \/l-e 2 



2GS sin /[(1+e cos /) cos 2 u— (1+e cos /o) cos 2 uo] 
c 2 na :i (l— e 2 ) 3 / 2 



GS{2[(/-/o)+e(sin/-sin/o)]-[(l+ecos/)sin2M-(l+ecos/o)sin2M ]} 
c 2 na 3 (l-e 2 ) 3 / 2 ' 



2GS cos j[2e(/-/ )+(l+e 2 )(sin /-sin /p)] 
c 2 nea 3 (l— e 2 ) 3 / 2 



2GS cos /(sin /—sin /o) 

r 2 r)pn 3 



(4) 

In eq. (jlj) n = ^jGM/a 3 is the unperturbed Keplerian mean motion, / is the 
true anomaly reckoning the instantaneous position of the test particle along 
the Keplerian ellipse, and u = uj + / is the argument of the latitude; it is 
intended that u = f + u since u is fixed for an unperturbed Keplerian orbit, 
where /o is the true anomaly at the epoch. Then, by means of the general 



Here R, T, N are the radial, transverse and out-of-plane (denoted also as normal) 
directions of the orthonormal frame co-moving with the test particle. Such a frame 
is usually adopted to decompose any perturbing acceleration acting on the particle 
itself. See eq. (JTSJ-eq. dJ) below. 



7 Similar results can be found in ISoffell (|l989l . p. 95), but they refer to the case 
/o = 0. 
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relations ( jCasotto 



19931 ). 



AR = (£) Aa-a cos /Ae + ^2s£AM, 
< AT = a sin / [l + ^2y] Ae + r(cos IAQ + Aw) + v / T 3 e2AA^, 



AN = r (sin ■uAJ — cos it sin I AQ) , 

(5) 

it is possible to analytically work out the Lense-Thirring perturbations of the 
R — T — N components of the orbit. They turn out to be 



AR LT = 2GScosI l 1 - cos (f-fo)] 



AT, 



LT 



VP 



2G5cosJ{2[(/-/o)-sin(/-/o)]+e[l-cos(/-/ )] sin/} 
c 2 na?\/l— e 2 (l+e cos /) 



(6) 



A at _ 2GSsin J{(l+ecos/o) cos up sin(/-/o)-[(/-/ )+e(sin /-sin/p)] cosu} 
LT c 2 na 2 VT^(l+e cos /) 



Note that the result of eq. ([6]), which is novel, is also an exact one; no approx- 
imations in e were used. Moreover, eq. ([6]) does not present any singularities 
for particular values of e and /. An inspection of eq. ([6]) shows that the radial 
perturbation ARi^ consists of the sum of a constant offset and a 1-cycle- 
per-revolution (cpr) harmonic term; no cumulative, secular components are 
present, so that the average radial shift is simply given by the constant term. 
Instead, in the transverse shift ATlt a dominant secular term is present in 
addition to a 1-cpr harmonic one; both are of zero order in e. There is also a 
smaller, harmonic component of order 0(e). The normal perturbation AA^lt 
has, at zero order in e, a secular term, whose amplitude is modulated by cosw, 
and a 1-cpr harmonic component; a smaller, harmonic term of order 0(e) is 
present as well. While ARuy and ATlt vanish for polar orbits, i.e. for / = 90 
deg, AA^lt is zero for equatorial orbits, i.e. for / = deg. 

In order to conveniently plot eq. (Hp as a function of time, we will use the 
useful relation for the guiding center (IMurray and Dermottl . Il999l . p. 41) 



/ ~ M + 2e sin M + \e 2 sin 2.M + 



+e 3 (g sin 3M - \ sin M ) + e 4 sin AM - g sin 2M ] 
indeed, M = n(t — t p ), where t p is the time of the passage at pericenter. 



(7) 
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The results of eq. (jSJ), with eq. (J7|), are the basic elements to work out 
the Lense-Thi rring perturbation Apu of the two-body range p. Indeed, from 
flCheneJ . liooih 

,2 



P 



p ' 

to be evaluated onto the unperturbed Keplerian ellipses 



flj(l-ej) 
1 + e j cos fj 



J 



A, B 



(9) 



of the two test particles A and B, it follows, for a generic perturbation (ICheng 

Ap = (Ar A - Ar B ) ■ p. 



20021 ). 



(10) 



In it 



Ar\ 



, - A/.', /? v + A/} Tj- + A.V, iVj, j = A, B, (11) 

where Ft, T, N are the unit vectors along the radial, transver se and out-of - 
plane directions, which are, in a body-centered {x, y, z} frame, (jChengi . 120021 ) 



R 



cos fl cos u — cos I sin Vt sin u 
sin Q cos it + cos I cos sin u 
sin / sin u 

- sin tt cos Q — cos / sin Q cos m 

- sin Q sin w + cos / cos Q cos w 

sin / cos u 



sin 1 sin i I 



(13) 



V 



7V — — sin I cos Q 
cos/. 

In the case of the unperturbed Keplerian ellipse it is 

r = r R, 



(14) 



/ 



(15) 



with r as in eq. and R given by eq. f[T2~j) . Concerning the gravitomagnetic 
field of the Earth, eq. (jSJ), evaluated for the spacecrafts A and B, has to be 
inserted into eq. (fT0|) -eq. (TTTj) to yield Ap LT ; eq. (JTj) allows to plot it as a 
function of time t. To avoid possible misunderstandings, it is important to 
point out that, in order to have the correct expression for the perturbation 
Ap L T of the intersatellite range, all the three components Ai? LT , AT LT , AiV LT 
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of eq. for both the satellites A and B concur to form Ap through eq. (ITU]) 
and eq. ( ITT]) . In other words, it would be incorrect to only consider Ai? L T and 
taking something like, say, Ap = ARa — AR B , although, at a first sight, one 
may be tempted to do so. 

Although more cumbersome, it is possible to analytically work out the two- 
body range-rate perturbation Ap (jCheng . 2002 ) as well. The first step consists 
of workin g out the R — T — N shifts of the test-particle's velocity. In general, 
they are ( jCasottol . Il993[ ) 



Av R = + ^) - na2 ^ = ^ (cos/Afi + Aco) - ^AM, 



Avt 



Av N 



na\J 1— e 2 
2r 



Aa + + ( cos /Afi + Aw ) > 



(l_ e 2)3/2 



[(cos u + e cos oS) AI + (sin u + e sin u) sin I Aft] 



(16) 



In the case of the gravitomagnetic field, the resulting velocity shifts are 



a LT _ 2GScosJ(l+ecos/)[2(/-/ )-sin(/-/o)+e(sin/-sin/ )] 



c 2 a 2 (l-e 2 ) 2 



. LT G5cos/{2(e+cos/)cos/ -(2+e 2 )+2(l+e 2 ) sin / sin f -e[2 cos /+4(/-/o) sin f-e cos 2/] } 
AV T = e 2 a 2 (l-e 2 ) 2 ' 



Av^ T = ^^^2)2 {2 (e cos bJ + cos m) [(1 + e cos / ) cos 2 u — (1 + e cos /) cos 2 



+ (e sin u; + sin u) [2 (/ — / ) + 2e (sin / — sin / ) — (1 + e cos /) sin 2u+ 



+ sin 2u + e cos f sin 2u ] } . 

(17) 

Also eq. ( ITT]) is an exact result in the sense that no approximations in e were 
used. 



Then, the follo wing unit vect or, computed onto the unperturbed Keplerian 
ellipse, is needed (jChengi . |2002| ) 



where 



{v A - vb) - P P 

J 

p 



P = (va - v B ) ■ p. 



(18) 
(19) 



It is intended that the Keplerian test particle's velocity has to be used in eq. 
(CBD-eq. (d; it is 

v = v R R + v T f, (20) 

with 

nae sin / 



vt 



VT- 



na(l+e cos /) 



(21) 



Note that, by constructio n, p„ is ortho gonal to p. The two-body range-rate 
perturbation Ap is, thus, (jChengl . 120021 1 



Ap = (Av A - Av B ) • p + (Ar A - Ar B ) • p u , 



where 



Aw, 



Av R Rj + Av J T Ti + Av 



Nn 



A, B. 



(22) 



(23) 



In the present case, inserting eq. ([6]) and eq. (fT7|) with eq. (J7J), allows to 
plot the gravitomagnetic range-rate perturbation as a function of time. Also 
in this case, we warn the reader that it would be incorrect to make, say, 
Ap = Av R — Av®: all the terms of eq. (jXTJ) for both satellites A and B must 
be used in eq. ( 1221 through eq. (l23il . An alternative approach to compute Ap 
consists of straightforwardly taking the derivative of Ap with respect to t after 
that its time series has been produced: indeed, it can be shown (see Figure [1] 
in Section I2TB1 that the result is the same. 



2.3 Discussion of the obtained results 



In FiguredJwe display the result of our numerical and analytical calculations 
for the gravitomagnetic field of the Earth; Table [2] resumes some quantitative 
features of the Lense-Thirring range and range-rate signatures. Concerning the 
analytical calculations for the range-rate displayed in the right-hand column of 
Figure [U the red signal (I) has been obtained by differencing the corresponding 
red range signal on the left. Instead, the green signal (II) has been obtained 
from the full application of the formalism dey eloped in Section 12.21 based on 
eq. (1T7|) . eq. ( 12"2"|) and eq. (12"31) ( jChengl . 120021 ) . As previously anticipated in 
Section [272| both the approaches are equivalent. Moreover, they agree with the 
numerical range-rate signal as well. At a first sight, one might find the pattern 
of the range-rate time series on the right somewhat puzzling if compared to 
that of the range ones on the left: there are intervals in the left-hand column 
where the plotted function is almost a straight line, while in the right-hand 
column its derivative is varying as rapidly as in other intervals. Actually, 
there is no contradiction at all, if one carefully looks at the scales displayed 
on the axes of the plots. Indeed, the magnitude of the range-rate signal is 
of the order of a few nm s -1 , which is perfectly in agreement with the quite 



10 



Lense-Thirring SST range perturbation: analytical calculation 



Lense-Thirring SST range-rate perturbation: analytical calculation (I) 




0.6 0.8 



Lense-Thirring SST range perturbation: numerical integration 
80 




Lense-Thirring SST range-rate perturbation: analytical calculation (II) 
6 
4 

" 2 
| 

5-2 

-4 

-6 




0.2 0.4 0.6 0.8 

1(d) 



0.4 0.6 0.8 

1(d) 




Fig. 1. First row from the top (red lines): analytically computed satellite-to-satellite 
range (left column) and range-rate (right column) signatures for GRACE A/B due 
to the Lense-Thirring effect. The range-rate signal has been obtained by differencing 
the range time series. Second row from the top (green line): analytically computed 
satellite-to-satellite range-rate (right column) signature for GRACE A/B due to 
the L ense-Thirring effect according to eq. (f22|) and eq. ([23]) with eq. (fT7|) (jCheng . 
2002). Third row from the top (blue lines): differences of the numerically integrated 
satellite-to-satellite range (left column) and range-rate (right column) signatures 
for GRACE A/B with and without the gravitomagnetic dynamical perturbation. 
The initial conditions, quoted in Table [IJ are common to both the perturbed and 
unperturbed integrations. The time span is AP = 1 d. The units are /j,m (range) 
and nm s" 1 (range-rate). 



small variations superimposed to the linear trend of the range, plotted with 
AP = 1 d = 86400 s. Moreover, from a visual inspection of Figure [1] it can 
be noticed that the range shift amounts to 80 /zm over 86400 s, yielding, on 
average, a range-rate of ~ 0.9 nm s -1 . Now, looking at the pictures in the 
right-hand column of Figure [IJ it appears that they are not exactly symmetric 
with respect to the horizontal axis, with an average slightly above it. This 
is fully confirmed by a quantitative numerical analysis of the range-rate time 
series of Figure [H which yields just an average of 0.9 nm s _1 . To further dispel 
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possible doubts on the reliability of our calculations, in Figure [2] we plot the 
analytical range and range-rate time series over a shorter time span, amounting 
to AP = (p| Kep) + Pjf ep) ) /2 = 0.065 d = 5630 s, so that it is easier to see 
how the changes in Ap are reflected by the time series of Ap. 

Lense-Thirring SST range perturbation: analytical calculation 

5 
4 

5 3 
<l 

2 
1 


0.01 0.02 0.03 0.04 0.05 0.06 
t(d) 

Lense-Thirring SST range-rate perturbation: analytical calculation (I) 
1.75 

1.5 

1.25 

7 „ 1 
g 

£ 0.75 
< 0.5 
0.25 


0.01 0.02 0.03 0.04 0.05 0.06 
t(ii) 

Fig. 2. First panel from the top: analytically computed satellite-to-satellite range 
signature, in pm, for GRACE A/B due to the Lense-Thirring effect. Second panel 
from the top: analytically computed satellite-to-satellite range-rate signature, in nm 
s , for GRACE A/B due to the Lense-Thirring effect. The details of the calculation 
are as in Figure [U The time span of the plots is AP = 0.065 d. 

As far as Figured] is concerned, the magnitude of the Lense-Thirring range 
signal is Aplt = 80 pm, while the gravitomagnetic range-rate effect is quite 
small, being of the order of Ap L T = 0.012 pm s" 1 (peak-to-peak amplitude). 
By assuming a present-day accuracy of a p ~ 5 pm in measuring the satellite- 
to-satellite range, the Lense-Thirring signature would fall, in principle, within 
the measurability domain at a ~ 6% level. Concerning the range-rate, if we 
assume a p ~ 0.5 pm s _1 , the Lense-Thirring effect on the range-rate Aplt is, 
instead, certainly too small to be detectable. 

Moving to the largest general relativistic orbital perturbation, i.e. the one 
due the gravitoelectric, Schwarzschild-like part of the gravitational field, Fig- 
ure [3] depicts its numerically integrated range and range-rate signatures. Table 
[2] summarizes some quantitative features of the gravitoelectric perturbations. 
The peak-to-peak amplitude of the range effect is as large as Aps crlw = 6000 
pm, so that it would be, in principle, detectable at a 8 x 10~ 4 level. The range- 
rate signal amounts to about Apg chw = 10 pm s _1 , which would be measurable 
with a nominal accuracy of about 5%. 
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Schwarzschild SST range perturbation: numerical integration Schwarzschild SST range-rate perturbation: numerical integration 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

t W) t (d) 

Fig. 3. Differences of the numerically integrated satellite-to-satellite range (left col- 
umn) and range-rate (right column) signatures for GRACE A/B with and without 
the Schwarzschild dynamical perturbation. The initial conditions, quoted in Table 
[U are common to both the perturbed and unperturbed integrations. The time span 
is AP = 1 d. The units are pm (range) and /jm s _1 (range-rate). 

Thus, GTR affects the orbital part of the GRACE satellite-to-satellite range 
in a detectable way, at least in principle, given the present-day level of accu- 
racy in measuring it. The same also holds for the range-rate, although only 
for the Schwarzschild signal, and at a lower level of accuracy with respect to 
the range. However, we wish to point out that it should, actually, be checked 

Table 2 

Peak-to-peak amplitudes of the inter-satellite GRACE range and range-rate pertur- 
bations Ap, A/9 caused by the general relativistic Schwarzschild and Lense-Thirring 
components of the Earth's gravitational field. The present-day accuracy in mea- 
suring the satellite-to satellite GRACE range and range-rate is a p < 10 pm and 
Op < 1 /jm s _1 , respectively. 



Dynamical effect 


Ap H 


Ap (/jm s 1 ) 


Schwarzschild 


6000 


10 


Lense-Thirring 


SO 


0.012 



if the relativistic signatures are not absorbed and removed from the range 
signal in estimating some of the various range parameters which are solved-for 
in the usual GRACE data processing. Indeed, it is exposed to some mismod- 
eled device behavior, which requires es timating many empir ical parameters in 
semi-dynamical orbit processing mode (jBalmino et al.l . 120051 ) . In fact, it would 
be necessary to realistically simulate the range and range-rate observations by 
fully modeling GTR, and implementing a data reduction: it is beyond the 
scope of the present work. On the other hand, our choice for the time inter- 
val of the numer ical integrations , close to the ones actually employed in real 
data reductions (jFlechtner et al.l . 1201 ll ). makes our analysis closer to what is 
actually done in real GRACE data analyses. 



In Section [3] we will look at the competing classical effects induced on 
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the range and range-rate of GRACE by some low-degree zonal harmonics 
J it & > 2 of the multipolar expansion of the Newtonian part of the Earth's 
gravitational potential accounting for its departures from spherical symmetry. 
Indeed, their unavoidably imperfect knowledge causes mismodeled range and 
range-rate signals which would corrupt the recovery of the relativistic ones at 
a level which has to be quantitatively assessed. On the other hand, such an 
investigation will also contribute to yield quantitative, although preliminary, 
evaluations of the level of a possible a-priori "imprinting" of GTR itself, not 
solved-for so far in all the GRACE-based models, on the estimated values of 
such zonals. This issue was treated, in the fram ework of the LAGEOS-based 
tests of the Lense-Thirring effect, by lloriol (120101 ) as far as the node precessions 
of the orbital planes of the GRACE spacecrafts are concerned. 



3 A-priori "imprint" level of GTR in the zonal KBR signature 



In this Section we numerically work out the effects on the GRACE range 
and range-rate caused by the mismodelling in the first six zonal coefficients 
of the terrestrial gravitational field. See Table [3] for their values and formal, 
statistical 1-a errors in one of the most recent global Earth's gravity field 
solution. 



Table 3 

Estimated values C^o and formal, statistical errors er^ of the normalized Stokes 
coefficients of the ge opotential f or i = 2, 3, 4, 5, 6, 7 from the GOCE/GRACE-based 
solution GOCO01S (|Pail et al.L [joifl l . Recall that J t = -^21 + 1 C l . 



Degree t 




°r 

w,o 




2 


-4.841649689 x 10" 4 


4 x 10" 


13 


3 


9.571980 x 10~ 7 


3 x 10" 


13 


4 


5.4000331 x 10" 7 


8 x 10" 


14 


5 


6.867018 x 10~ 8 


6 x 10" 


11 


G 


-1.4995817 x 10" 7 


4 x 10" 


14 


7 


9.051062 x 10" 8 


4 x 10" 


14 
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3.1 The even zonals 



The range and range-rate perturbations due to the mismodelling in the first 
three even zonals are depicted in Figure H] and quantitatively summarized in 
Table H 



J 2 SST range perturbation: numerical integration 
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£ -40 
-50 
-60 




0.2 0.4 0.6 

1(d) 



J 4 SST range perturbation: numerical integration 










-5 




-10 
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-15 
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J 6 SST range perturbation: numerical integration 




J 2 SST range-rate perturbation: numerical integration 




0.4 0.6 
1(d) 



J 4 SST range-rate perturbation: numerical integration 




0.4 0.6 
1(d) 



J 6 SST range-rate perturbation: numerical integration 




Fig. 4. Differences of the numerically integrated ranges (left column) and range-rates 
(right column) for GRACE A/B with and without the classical dynamical perturba- 
tions due to the mismodelled even zonal harmonics Ji = —\JU + 1 Ce o, I = 2, 4, 6 . 



Acco rding to the recent combined GOCE-GRACE solution GOCO01S (jPail et al. 
the formal, statistical uncertainties in the normalized Stokes coefficients are, 

a-r, = 0.44 x lCT 12 ,^ = 0.8 x 1CT 13 , 077 = 0.4 x 10" 13 respectively. The initial 

02,0 ' c.4,0 ' c-6,0 1 J 

conditions, quoted in Table [U are common to both the perturbed and unperturbed 
integrations. The time span is AP = 1 d. The units are /im. (range) and /mi s _1 
(range-rate) . 



According to Table |2] and Table HI the mismodelled ranges due to the even 
zonals are slightly smaller than the corresponding Lense-Thirring effect, while 
they are about 90 — 600 times smaller than the Schwarzschild range; instead, 
the Schwarzschild range-rate signal is 700 — 4000 times larger than the cor- 
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Table 4 

Peak-to-peak amplitudes of the inter-satellite GRACE range and range-rate pertur- 
bations Ap, Ap caused by the first three mismodelled even zonal harmonics of the 
classical part of the geopotential over AP = 1 d. The figures for , = 2,4, 6 
along with their nominal values, can be found in Table [3j 



Dynamical effect 


Ap H 


Ap (pm s x ) 


■h 


70 


0.014 


Ji 


30 


0.006 


J 6 


10 


0.0025 



responding mismodelled signatures due to the even zonals. However, it must 
be considered that, given a specific Earth's gravity field model, the actual un- 
certainties in its estimated even zonals may be up to one order of magnitude 
larger with respect to their formal, statistical errors. Moreover, an even more 
conservative approach to realistically evaluate the true uncertainties in the 
even zonals consists of comparing their estimated v alues from different global 



gravity field solutions fjWagner and McAdool . 120121 ) 



Conversely, we can use Table [2] and Table H] (or, equivalently, Figure [Ql3] and 
Figure H|) to obtain preliminary, quantitative evaluations of a possible a-priori 
"imprinting" of GTR itself in the even zonals considered. By posing x = p, p, 
it can be done from 



Ax GTR = A&A , = (24) 

Jo 



Thus, by dividing the relativistic range and range-rate perturbations Ax GTR 
by the corresponding normalized classical ones A£g for each degree I considered 
gives us a sort of "effective" relativistic even zonal Ajj; eS \ i.e. the part of the 
even zonal of degree £ which would give a signal as large as those due to GTR. 

According to Table [5] and Table and in view of the present-day level of 
accuracy in estimating them (see Table[3]), the a-priori "imprinting" of GTR on 
the even zonals should not be neglected, especi ally a s far as the Schwarzschild 



part is concerned. Following the approach by lloriol (120101 ). it can be shown 
that a Lense-Thirring "imprint" in J 4 and Jq as large as that of Table [5]- Table 
[6] would correspond to a rs 0.2% systematic bias in the LAGEOS-based tests 
of such a relativistic effect. However, there is room for a "contamination" of 
GTR itself due to the gravitoelectric part. 
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Table 5 

"Effective" general relativistic parts Ajj, c ^ of the even zonals of degree I = 2, 4, 6 
obtained from the range perturbations Ap. They are a preliminary measure of a pos- 
sible a-priori "imprinting" of GTR itself, not explicitly solved-for in all the GRACE- 
based solutions produced so far, on the even zonals. Compare these figures with the 
formal, statistical errors <j-tj , t = 2,4, 6 in Table O 



Dynamical effect (Ap) 


Aj( eff) 


Aji cff) 


a4 cS) 


Schwarzschild 


3 x KT 11 


2 x KT 11 


2 x KT 11 


Lense-Thirring 


5 x KT 13 


2 x KT 13 


3 x KT 13 



Table 6 

"Effective" general relativistic parts A Jg°^ of the even zonals of degree i = 2,4,6 
obtained from the range-rate perturbations Ap. They are another preliminary mea- 
sure of a possible a-priori "imprinting" of GTR itself, not explicitly solved-for in 
all the GRACE-based solutions produced so far, on the even zonals. Compare these 
figures with the formal, statistical errors o-q , i = 2,4, 6 in Table El 



Dynamical effect (Ap) 


a4 cS) 


Ajf^ 


AJ^ 


Schwarzschild 
Lense-Thirring 


3 x KT 10 
3 x KT 13 


1 x KT 10 

2 x KT 13 


2 x KT 10 
2 x KT 13 



3.2 The odd zonals 

In Figure [5] and Table [7] we repeat the same analysis for the first three odd 
zonals. 

Table 7 

Peak-to-peak amplitudes of the inter-satellite GRACE range and range-rate pertur- 
bations Ap, Ap caused by the first three mismodelled odd zonal harmonics of the 
classical part of the geopotential over AP = 1 d. The figures for <7q , £ = 3,5,7 
along with their nominal values, can be found in Table El 



Dynamical effect 


Ap (pm) 


Ap (pm s x ) 


■h 


70 


0.015 


■h 


20 


0.004 


■h 


0.5 


0.0004 
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J3 SST range perturbation: numerical integration 



J 3 SST range-rate perturbation: numerical integration 
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J5 SST range perturbation: numerical integration 
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J 7 SST range perturbation: numerical integration 
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J5 SST range-rate perturbation: numerical integration 
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J 7 SST range-rate perturbation: numerical integration 
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Fig. 5. Differences of the numerically integrated ranges (left column) and range-rates 
(right column) for GRACE A/B with and without the classical dynamical perturba- 
tions due to the mismodelled odd zonal harmonics Ji = —\J1l + 1 C^o, I = 3, 5, 7 . 



Acco rding to the recent combined GOCE-GRACE solution GOCO01S (jPail et al. 



the formal, statistical uncertainties in the normalized Stokes coefficients are, 
_ n o7 w 1 n-12 „_ _ 0.6 x 1Q- 13 , o-^ 7 q = 0.4 x 10~ 13 respectively. The initial 



oj, =0.27x10 

L-3,0 



L-5,0 



conditions, quoted in Tabled! are common to both the perturbed and unperturbed 
integrations. The time span is AP = 1 d. The units are /mi (range) and /mi s _1 
(range-rate) . 



Concerning the range, it turns out that the shifts due to the mismodelled 
odd zonals are about as large as the Lense-Thirring signature and smaller 
than the Schwarzschild one by a factor pa 85 — 12000. The gravitomagnetic 
range-rate perturbation is of the same order of magnitude of the J3— induced 
range-rate mismodelled signature, but it is larger than the other ones up to 160 
times (J7). The Schwarzschild range-rate effect is larger than the competing 
odd zonals ones by a factor w 670 — 25000. The Lense-Thirring range-rate 
signal is smaller than the mismodelled J3 one, but it is 3 — 30 times larger 
than the J5 and J-j mismodelled signatures. 

The a-priori "imprinting" of GTR on the odd zonals can preliminarily be 
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quantified according to Table [8] (range) and Table [9] (range-rate). 
Table 8 

"Effective" general relativistic parts A jj cff ^ of the odd zonals of degree t = 3,5,7 
obtained from the range perturbations Ap. They are a preliminary measure of a pos- 
sible a-priori "imprinting" of GTR itself, not explicitly solved-for in all the GRACE- 
based solutions produced so far, on the odd zonals. Compare these figures with the 
formal, statistical errors <jtj , I = 3, 5, 7 in Table [3J 



Dynamical effect (Ap) 


a4 cS) 


a4 cS) 




Schwarzschild 
Lense-Thirring 


3 x KT 11 

4 x KT 13 


2 x KT 11 
2 x 1CT 13 


5 x KT 10 
6.4 x KT 12 



Table 9 

"Effective" general relativistic parts Ajj eff ^ of the odd zonals of degree £ = 3,5,7 
obtained from the range-rate perturbations Ap. They are another preliminary mea- 
sure of a possible a-priori "imprinting" of GTR itself, not explicitly solved-for in 
all the GRACE-based solutions produced so far, on the odd zonals. Compare these 
figures with the formal, statistical errors o-q , i = 3, 5, 7 in Tabled 



Dynamical effect (Ap) 


Aj( cff ) 


a4 cS) 


A4 cff ) 


Schwarzschild 
Lense-Thirring 


2 x KT 10 
2 x KT 13 


1 x KT 10 

2 x KT 13 


1 x KT 9 
1 x KT 12 



Also in this case, the present-day level of accuracy in determining the odd 
zonal coefficients of the geopotential does not allow, in principle, to neglect 
GTR SIS db potential source of a-priori aliasing. 

As far as the LAGEOS-based tests of the Lense-Thirring effect are con- 
cerned, the issue of the a-priori "imprinting" of GTR on the odd zonals is not 
a concern since the node of a satellite is not secularly affected by them. 



4 The impact of the time delay on the GRACE SST range 

In this Section we treat the range shifts Ap coming from the GTR-induced 
time delays At affecting the intersatellite range. 
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4-1 The gravitomagnetic time delay 



Here we deal in detail with the effect of the gravitomagnetic field of the 
Earth on the propagation of the electromagnetic waves between both GRACE 
A and GRACE B, which is not modeled in the softwares used to process 
GRACE data. 



The Len se-Thirring tim e delay A^lt between A and B can conveniently be 
written as (jCiufolini et al.l . 120031 ) 



At 



2GS 



LT 



1 1 



S • (f A X f B ) 
1 + r A • f B 



(25) 



Generally speaking, for a pair of satellites orbiting at the altitude of GRACE 
it is 



2GS 



1 1 



2.8 x 1(T 17 s, 



(26) 



corresponding to a range shift of just 8.5 nm. Moreover, in the particular case 
of GRACE, the geometric factor entering eq. ( 1251) is quite small because both 
r A and r B almost lie in a plane containing S as well. Indeed, it can be showed 
that 

S • (f A X r B ) 



~ HT' 3 - 10" 



(27) 



Thus, we conclude that the part of Ap due to the Lense-Thirring time delay 
is completely negligible for the GRACE mission. 



4-2 The Shapiro time delay 



Concerning the impact of the gravitoelectric part of the terrestrial gravita- 
tional field on the propagation of the electromagnetic waves linking GRACE 
A and GRACE B, and included in the GRACE models, the Shapiro time delay 
fjShapirol . E96l affect ing the SST range can be valuated as[f] (jPetit and Luzuml . 
2010|, p. 164) 



At 



2GM 



GE 



In 



ta + t b + p 



(28) 



r A + r B -pJ 

It can be shown that eq. (128]) corresponds to a secularly increasing range shift 
cA^ge ~ 200 pm over AP = 1 d. Thus, although the Shapiro delay on the 



See also (jMoverl . 120031 ). 
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SSR range is detectable, it is about 30 times smaller than the range shift due 
to the orbital motions. 



5 Summary and conclusions 

Given the present-day high level of accuracy in measuring the GRACE 
satellite-to-satellite biased range [p p < 10 /jm) and range-rate [p p < 1 //m 
s -1 ), we preliminarily investigated the impact of GTR on such directly observ- 
able quantities. We did not consider the general relativistic effects connected 
with the propagation of the electromagnetic waves linking the two spacecrafts. 
A further motivation for the present sensitivity analysis is given by the cur- 
rently ongoing efforts to design a follow-on of GRACE accurate to nm s -1 , or 
better, in the range-rate. Moreover, GTR has never been solved-for in all the 
GRACE-based Earth's global gravity field solutions produced so far, so that 
the multipoles of the terrestrial gravitational field estimated in them may, in 
principle, retain an a-priori "imprinting" of GTR itself. This fact is important 
since some of such Earth's global gravity models were used as background 
reference models in some tests of GTR itself performed with other satellites. 

By numerically integrating the GRACE A/B post-Newtonian equations of 
motion in a geocentric frame over a time span 1 d long, corresponding to about 
15 full orbits, we found that the GTR range signals are as large as 6000 /im 
(Schwarzschild) and 80 [xm (Lense-Thirring), while the sizes of the range-rate 
GTR effects are 10 /xm s" 1 (Schwarzschild) and 0.012 fim s _1 (Lense-Thirring). 
An analytical calculation for the Lense-Thirring effect confirmed the numer- 
ical results concerning it. If, on the one hand, they are larger than <j p and 
(jp, apart from the Lense-Thirring range-rate, on the other hand the imperfect 
knowledge of some low-degree zonal harmonics of the geopotential causes com- 
peting range and range-rate perturbations which would corrupt the recovery 
of the relativistic signals of interest. According to the formal, statistical errors 
released in one of the latest GOCE/GRACE-based models, the mismodelled 
signatures of the first six zonals are smaller than the Schwarzschild range and 
range-rate ones by several orders of magnitude. Instead, they are about of the 
same order of magnitude of, or slightly smaller than, the Lense-Thirring range 
and range-rate perturbations. Conversely, by comparing the relativistic and 
the zonal orbital effects on the GRACE range and range-rate it was possible to 
quantitatively assess the level of a possible a-priori "imprinting" of GTR itself 
on the zonals considered. It turned out that it is not negligible as far as the 
Schwarzschild component is concerned, while the Lense-Thirring "imprint" is 
about at the edge of the present-day level of accuracy in determining them. 

As a caveat concerning the actual measurability of the GTR effects consid- 
ered here, we stress the need of checking with extensive numerical simulations 
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if the relativistic signatures are not absorbed and removed from the range 
signal in estimating some of the various range parameters which are solved-for 
in the usual GRACE data processing. Implementing such a non-trivial task is 
outside the scopes of the present paper, but it could-and, in fact, should-be 
done by the various groups worldwide engaged in analyzing longer and longer 
GRACE data records. It is highly desirable that they produce dedicated global 
gravity field models explicitly estimating GTR as well along with the usual 
Stokes coefficients of the geopotential. 

Finally, let us note that the approach presented here can, in principle, 
also be extended to other satellite-to-satellite orbital configurations suitably 
designed to enhance the relativistic signatures, and to the dynamical effects 
caused by various modified models of gravity. Such a possibility has recently 
become appealing after the proposal of a putative Earth-sourced fifth force as 
a possible explanation of Earth-based neutrino phenomenology. 
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